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Abstract. A simple model accounting for the ejection of heavy particles from the 
vortical structures of a turbulent flow is introduced. This model involves a space and 
time discretization of the dynamics and depends on only two parameters: the fraction 
of space-time occupied by rotating structures of the carrier flow and the rate at which 
particles are ejected from them. The latter can be heuristically related to the response 
time of the particles and hence measure their inertia. It is shown that such a model 
reproduces qualitatively most aspects of the spatial distribution of heavy particles 
transported by realistic flows. In particular the probability density function of the 
mass m in a cell displays an power-law behavior at small values and decreases faster 
than exponentially at large values. The dependence of the exponent of the first tail 
upon the parameters of the dynamics is explicitly derived for the model. The right tail 
is shown to decrease as exp(— Cm log m). Finally, the distribution of mass averaged 
over several cells is shown to obey rescaling properties as a function of the coarse-grain 
size and of the ejection rate of the particles. Contrarily to what has been observed 
in direct numerical simulations of turbulent flows (Bee et at, nlin. CD/0608045), such 
rescaling properties are only due in the model to the mass dynamics of the particles 
and do not involve any scaling properties in the spatial structure of the carrier flow. 



PACS numbers: 47.55.Kf, 05.50.+q, 47.27.T- 
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1. Introduction 

Understanding the dynamics of small-size tracer particles or of a passive field transported 
by an incompressible turbulent flow plays an important role in the description of several 
natural and industrial phenomena. For instance it is well known that turbulence has 
the property to induce an efficient mixing over the whole range of length and time 
scales spanned by the turbulent cascade of kinetic energy (see e.g. PQ). Describing 
quantitatively such a mixing has consequences in the design of engines, in the prediction 
of pollutant dispersion or in the development of climate models accounting for transport 
of salinity and temperature by large-scale ocean streams. 

However, in some settings, the suspended particles have a finite size and a mass 
density very different from that of the fluid. Thus they can hardly be modeled by 
tracers because they have inertia. In order to fully describe the dynamics of such 
inertial particles, one has to consider many forces that are exerted by the fluid even 
in the simple approximation where the particle is a hard sphere much smaller than the 
smallest active scale of the fluid flow [2] . Nevertheless the dynamics drastically simplifies 
in the asymptotics of particles much heavier than the carrier fluid. In that case, and 
when buoyancy is neglected, they interact with the flow only through a viscous drag, so 
that their trajectories are solutions to the Newton equation : 



d 2 X 1 



alt 2 t 
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where u denotes the underlying fluid velocity field and r is the response time of the 
particles. Even if the carrier fluid is incompressible, the dynamics of such heavy particles 
lags behind that of the fluid and is not volume-preserving. At large times particles 
concentrate on singular sets evolving with the fluid motion, leading to the appearance of 
strong spatial inhomogeneities dubbed preferential concentrations. At the experimental 
level such inhomogeneities have been known for a long time (see [3] for a review). At 
present the statistical description of particle concentration is a largely open question 
with many applications. We mention the formation of rain droplets in warm clouds [I], 
the coexistence of plankton species [5], the dispersion in the atmosphere of spores, dust, 
pollen, and chemicals [6], and the formation of planets by accretion of dust in gaseous 
circumstellar disks [7]. 

The dynamics of inertial particles in turbulent flows involves a competition between 
two effects: on the one hand particles have a dissipative dynamics, leading to the 
convergence of their trajectories onto a dynamical attractor [8], and on the other hand, 
the particles are ejected from the coherent vortical structures of the flow by centrifugal 
inertial forces [9]. The simultaneous presence of these two mechanisms has so far led to 
the failure of all attempts made to obtain analytically the dynamical properties or the 
mass distribution of inertial particles. In order to circumvent such difficulties a simple 
idea is to tackle independently the two aspects by studying toy models, either for the 
fluid velocity field, or for the particle dynamics that are hopefully relevant in some 
asymptotics (small or large response times, large observation scales, etc.). Recently 
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an important effort has been made in the understanding of the dynamics of particles 
suspended in flows that are 5-correlated in time, as in the case of the well-known 
Kraichnan model for passive tracers (TO]. Such settings, which describe well the limit 
of large response time of the particles, allows one to obtain closed equations for density 
correlations by Markov techniques. The ^-correlation in time, of course, rules out the 
presence of any persistent structure in the flow; hence any observed concentrations can 
only stem from the dissipative dynamics. Most studies in such simplified flows dealt 
with the study of the separation between two particles pj], [12], [13j El 03] • 

Recent numerical studies in fully developped turbulent flows [16] showed that the 
spatial distribution of particles at lengthscales within the inertial range are strongly 
influenced by the presence of voids at all active scales spanned by the turbulent cascade 
of kinetic energy. The presence of these voids has a noticeable statistical signature on 
the probability density function (PDF) of the coarse-grained mass of particles which 
displays an algebraic tail at small values. To understand at least from a qualitative 
and phenomenological viewpoint such phenomena, it is clearly important to consider 
flows with persistent vortical structures which are ejecting heavy particles. For this 
purpose, we introduce in this paper a toy model where the vorticity field of the carrier 
flow is assumed piecewise constant in both time and space and takes either a finite fixed 
value u or vanishes. In addition to this crude simplification of the spatial structure of 
the fluid velocity field we assume that the particle mass dynamics obeys the following 
rule: during each time step there is a mass transfer between the cells having vorticity uj 
toward the neighboring cells where the vorticity vanishes. The amount of mass escaping 
to neighbors is at most a fixed fraction 7 of the mass initially contained in the ejecting 
cell. We show that such a simplified dynamics reproduces many aspects of the mass 
distribution of heavy particles in incompressible flow. In particular, we show that the 
PDF of the mass of inertial particles has an algebraic tail at small values and decreases 
as exp(— Am logm) when m is large. Analytical predictions are confirmed by numerical 
experiments in one and two dimensions. 

In section [2] we give some heuristic motivations for considering such a model and 
a qualitative relation between the ejection rate 7 and the response time r of the heavy 
particles. Section [3] consists in a precise definition of the model in one dimension and in 
its extension to higher dimensions. Section @] is devoted to the study in the statistical 
steady state of the PDF of the mass in a single cell. In section [5] we study the mass 
distribution averaged over several cells to gain some insight on the scaling properties in 
the mass distribution induced by the model. Section [6] encompasses concluding remarks 
and discussions on the extensions and improvements of the model that are required to 
get a more quantitative insight on the preferential concentration of heavy particles in 
turbulent flows. 
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2. Ejection of heavy particles from eddies 



The goal of this section is to give some phenomenological arguments explaining why the 
model which is shortly described above, might be of relevance to the dynamics of heavy 
particles suspended in incompressible flows. In particular we explain why a fraction of 
the mass of particles exits a rotating region and give a qualitative relation between the 
ejection rate 7 and the response time r entering the dynamics of heavy particles. For 
this we focus on the two-dimensional case and consider a cell of size £ where the fluid 
vorticity u is constant and the fluid velocity vanishes at the center of the cell. This 
amounts to considering that the fluid velocity is linear in the cell with a profile given 
to leading order by the strain matrix. Having a constant vorticity in a two-dimensional 
incompressible flow means that we focus on cases where the two eigenvalues of the strain 
matrix are purely imaginary complex conjugate. The particle dynamics reduces to the 
second-order two-dimensional linear system 
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It is easily checked that the four eigenvalues of the evolution matrix are the following 
complex conjugate 
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This means that the distance of the particles to the center of the cell increases 
exponentially fast in time with a rate \i. If we now consider that the particles are initially 
uniformly distributed inside the cell, we obtain that the mass of particles remaining in 
it decreases exponentially fast in time with a rate equal to — 2/i. Namely the mass of 
particles which are still in the cell at time T is 

T / 1 



m(T) = m(0) (1 — 7) = m(0) exp 



— -1 + 



2Vl + 16r 2 



.(5) 



The rate 7 at which particles are expelled from the cell depends upon the response time 
t of the particles and upon two characteristic times associated to the fluid velocity. The 
first is the time length T of the ejection process which is given by the typical life time of 
the structure with vorticity u. The second time scale is the turnover u^ 1 which measures 
the strength of the eddy. There are hence two dimensionless parameters entering the 
ejection rate 7: the Stokes number St = tuj giving a measure of inertia and the Kubo 
number Ku = Tuj which is the ratio between the correlation time of structures and their 
eddy turnover time. One hence obtain the following estimate of the ejection rate 
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The graph of the fraction of particles ejected from the cell as a function of the Stokes 
number is represented in figure [1] for three different values of the Kubo number. The 
function goes to zero as KuSt in the limit St — > and as KuSt~ l l 2 in the limit St — > oo. 
It reaches a maximum which is an indication of a maximal segregation of the particles, 
for St « 1.03 independently of the value of Ku. 

st=\m 



Ku = 3 




Stokes number St = x CO 



Figure 1. Fraction of the mass of particles that are uniformly distributed and ejected 
from an eddy of arbitrary size £ as a function of the Stokes number St = tlj. The 
various curves refer to different values of the Kubo number Ku as labeled. 

In three dimensions, one can extend the previous approach to obtain an ejection 
rate for cells with a uniform rotation, i.e. a constant vorticity uj. There are however 
two main difficulties. The first is that in three dimensions the eigenvalues of the strain 
matrix in rotating regions are not anymore purely imaginary but have a real part given 
by the opposite of the rate in the stretching direction. Such a vortex stretching has to 
be considered to match observation in real flows. The second difficulty stems from the 
fact that the vorticity is now a vector and has a direction, so that ejection from the 
cell can be done only in the directions perpendicular to the direction of uj. These two 
difficulties imply that the spectrum of possible ejection rates is much broader than in 
the two-dimensional case. However the rough qualitative picture is not changed. 

3. A simple mass transport model 

We here describe with details the model in one dimension and mention at the end of the 
section how to generalize it to two and higher dimensions. Let us consider a discrete 
partition of an interval in N small cells. Each of these cell is associated to a mass which 
is a continuous variable. We denote by rrij(n) the mass in the j-th cell at time t = n. At 
each integer time we choose randomly N independent variables; fL = 1 with probability 
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p and Qj = with probability 1 — p. The evolution of mass between times n and n + 1 
is given by: 



m. 



(n+1) 



rn 



in) [2 - £lj_x - Q j+ i] mj(n) 



if a 
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(7) 



In other terms, when fL = 1, the j-th cell looses mass if Q 



or fL 



0, and 

1. The flux of mass between the 
j-th and the (j + l)-th cell is proportional fL, — Qj+i (see figure [2]). In particular, if 



when fL, 



0, it gains mass if Qj-i = 1 or Q 



■3+1 



1 - 



Ci 



Figure 2. Sketch of the dynamics in the one-dimensional case: the fluxes of mass are 
represented as arrows. A cross means no flux. 

Qj = Qj+%, no mass is transfered between cells. When the system is supplemented by 
periodic boundary conditions between the cells N and 1, it is clear that the total mass is 
conserved. Hereafter we assume that the mass is initially nij = 1 in all cells., so that the 
total mass is ^ • rrij = N. Spatial homogeneity of the random process Qj implies that 
{rrij) = 1 for all later times, where the angular brackets denote average with respect to 
the realizations of the fL/s. 

A noticeable advantage of such a model for mass transportation is that the mass field 
m = (mi, . . . , m^) defines a Markov process. Its probability distribution p N (m, n + 1) 
at time n + 1, which is the joint PDF of the masses in all cells, is related to that at time 
n by a Markov equation, which under its general form can be written as 



p N (m,n + 1) 



al N m! ppf(m' ', n) P[m' — > m] 



J d N m'p N (m',n) J d N Q p(Sl) P[m' -> m\Q] 



(8) 



where P[m' — > m|f2] denotes the transition probability from the field m' to the field 
m conditioned on the realization of f2 = (Q\, . . . , Qn)- I n our case it takes the form 

N 

P[m' -> m\Q,] = Y[ 5[rrij - (m'j + ^-i(n) - Hj{n))\ . (9) 

3=1 

The variable p,j denotes here the flux of mass between the j-th and the (j + l)-th cell. 
It is a function of Qj, Qj+x, and of the mass contained in the two cells. It can be written 

as 

Hj{n) = | [Qj(n)(l - Qj +1 (n))m'j{n) - Q j+1 {n)(l - Qj(n))m' j+1 (n)] . (10) 

In the particular case we are considering, the joint probability of the fL/s factorizes and 
we have 



p(Qj) =p8{Qj - 1) + (1 -p) 8{Qj) 



ID 
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so that the Markov equation (jSJ) can be written in a rather explicit and simple manner. 

The extension of the model to two dimensions is straightforward. The mass transfer 
out from a rotating cell can occur to one, two, three or four of its direct nearest neighbors 
(see figure [3] left). One can similarly derive a Markov equation which is similar to (JH]) 
for the joint PDF pjv,iv(M, n) at time n of the mass configuration M = {^i,j}i<i,j<v- 
The transition probability reads in that case 

N N 



Ifil 



run 

t=i j=i 

where the fluxes now take the form 

7 
4 



P[P 



m 



m itj + fi 



(i) 



u (1) ^ 



+ A*S-i-/*S)].(12) 



(2)^ 



(13) 



and //^- defined by inverting i and j in the definition of . 

After a large number of time steps, a statistically steady state is reached. The 
stationary distribution is obtained assuming that p NtN (M, n) = p N)N (M) is independent 
of n in the Markov equation (jHJ). In this stationary state the mass fluctuates around its 
mean value 1 corresponding to a uniform distribution; strong deviations at small masses 
can be qualitatively observed (see figure [3] right). 
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Figure 3. Left: sketch of the ejection model in two dimensions; the rotating cells are 
designated by small eddies; the flux of mass (blue arrows) is from the rotating cells to 
those without any eddy. Right: snapshot of the distribution of mass in the statistically 
steady regime for a domain of 50 2 cells with p — 0.75; white squares are almost empty 
cells and the darker regions correspond to cells where the mass is concentrated. 



The model can be easily generalized to arbitrary space dimension. However, as 
we have seen in previous section, besides its interest from a purely theoretical point of 
view, the straightforward extension to the three-dimensional case might not be relevant 
to describe concentrations of inertial particles in turbulent flows. 



4. Distribution of mass 



Let us consider first the one-dimensional case in the statistically stationary regime. 
After integrating (jSJ), one can express the single-point mass PDF pi in terms of the 
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three-point mass distribution p 3 at time n 

PI K> = /i M ;.>»;i< lPj « 1 , m ;,,„; t ,)/« J „ 1 « J a J ,, * 

p{flj-i) p(flj) p(fy+i) 8[mj - (m'j + Hj^i - /ij)] . (14) 

We then explicit all possible fluxes, together with their probabilities, by considering 
all possible configurations of the spin vorticity triplet (fij_i, fl/, fij+i). The results are 
summarized in table [TJ This leads to rewrite the one-point PDF as 

Pi(m)= [p 3 + (l-p) 3 ] pi(m) + 2] j pi( j _™ /2 ) + 

, P(l-P) 2 / m \ , a [ 2 ™ h , ( , 7 A 

H : pi + 2p (1 - p) dm p2\m ,m - -m + 

1-7 V 1 - 7/ Jo v 2 / 

/■2m/7 p2m/-y-m' 

+ p 2 (l — p) / dm'/ dm'' p 3 ym',m — —(m' + m''),m'j . (15) 



'o v/o 

The first term on the right-hand side comes from realizations with no flux. The second 
term is ejection to one neighbor and the third to two neighboring cells. The fourth term 
involving an average over the two-cell distribution is related to events when mass is 
transfered from a single neighbor to the considered cell. Finally, the last term accounts 
for transfers from the two direct neighbors. Note that, in order to write ( [TBI) , we made 
use of the fact that P2(x, y) = P2^y, x). 



Table 1. Enumeration of all possible configurations of the spin vorticity fl in three 
neighboring cells, together with their probabilities and the associated mass fluxes. 
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Numerical simulations of this one- dimensional mass transport model are useful to 
grab qualitative information on p\. Figure H] represents the functional form of pi in 
the stationary regime for various values of the ejection rate 7 and for p = 1/2. The 
curves are surprisingly similar to measurements of the spatial distribution of heavy 
particles suspended in homogeneous isotropic turbulent flows [3l[T6]. This gives strong 
evidence that, on a qualitative level, the model we consider reproduces rather well the 
main mechanisms of preferential concentration. More specifically, a first observation is 
that in both settings the probability density functions display an algebraic behavior for 
small masses. This implies that the ejection from cells with vorticity one has a strong 
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Figure 4. Log- log plot of the one-point PDF of the mass in one dimension for p = 1/2 
and different values of the parameter 7 as labeled. The integration was done on a 
domain of 2 16 = 65536 cells and time average wee performed during 10 6 time steps 
after a statistical steady state is reached. 



statistical signature. A second observation is that at large masses, the PDF decays 
faster than exponentially, as also observed in realistic flows. As we will now see these 
two tails can be understood analytically for the model under consideration. 

We here first present an argument explaining why an algebraic tail is present at 
small masses. For this we exhibit a lower bound of the probability P < (m) that the mass 
in the given cell is less than m. Namely, we have 

P < (m) = Prob(ra,-(n) < m) > Prob (A) , (16) 

where A is a set of space-time realizations of Q such that the mass in the j-th cell at 
time n is smaller than m. For instance we can choose the set of realizations which are 
ejecting mass in the most efficient way: during a time N before n, the j-th cell has spin 
vorticity 1 and its two neighbors have 0. The mass at time n is related to the mass at 
time n — N by 

/ s xn , log[m,-(n — N)/mJn)] 

mj {n) = (1 - 7) A m j (n - N) , that is iV = 61 y — J j- n n . (17) 

The probability of such a realization is clearly ^(1 — p) 2N . Replacing iV by the 
expression obtained in (fT?I) . we see that 



Prob (A) 



rrij(n) 



i/ 3 



rrij (n — N) 



with (3 



log[p(l - pf 



log(l-7) 

After averaging with respect to the initial mass rrij(n — N), one finally obtains 

P < (m)>Am> 3 . (19) 
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Hence the cumulative probability of mass cannot have a tail faster than a power law at 
small arguments. It is thus reasonable to make the ansatz that pi(m) have an algebraic 
tail at m — > 0, i.e. that pi(m) ~ Cm a . To obtain how the exponent a behaves as a 
function of the parameters 7 and p, this ansatz is injected in the stationary version 
of the Markov equation (1151) . One expects that the small- mass behavior involves only 
the terms due to ejection from a cell, namely the three first terms in the r.h.s. of ffT5j) . 
and that the terms involving averages of the two-point and three-point PDFs give only 
sub-dominant contributions. This leads to 

2p 2 (l — p) f m 



Cm a 



C[p 3 + (l-p) 3 ] 



m 
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7/2 
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rn 
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Equating the various constants we finally obtain that the exponent a satisfies 
2p (1-p) 



+ 



3. 



(20) 



(21) 



(l-7/2) Q+1 (1-7) Q+1 

Note that the actual exponent a given by this relation is different from the lower-bound 
(3 + 1 obtained above in ffl8|) and f|T9|) . However it is easily checked that a approach 
the lower bound when p — > 0. As seen from figure EJ formula ( |2TI) is in good agreement 
with numerics. Note that the large error bars obtained for p small and 7 large are due 




0.4 0.6 

Ejection rate y 

Figure 5. Scaling exponent a as a function of the ejection rate 7 for three different 
values of p as labeled. The solid lines represents the prediction given by l[2T|) : the error 
bars are estimated from the maximal deviation of the logarithmic derivative from the 
estimated value. Inset: difference between the numerical estimation and the value 
predicted by theory. 



to the presence of logarithmic oscillations in the left tail of the PDF of mass. This log 
periodicity is slightly visible for 7 = 0.9 in figure HI It occurs when the spreading of the 
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distribution close to the mean value m = 1 is much smaller than the rate at which mass 
is ejected. This results in the presence of bumps in the PDF at values of m which are 
powers of 1 — 27. Notice that for all values of p, one has a < when 7 > 2/3. However, 
according to the estimate (|SJ), values of the ejection rate larger than 2/3 can be attained 
only for large enough Kubo numbers. This is consistent with the fact that power-law 
tails with a negative exponent were not observed in the direct numerical simulations of 
turbulent fluid flows [16] where Ku ~ 1. 

It is much less easy to get from numerics the behavior of the right tail of the mass 
PDF pi(m). As seen from figure HI there was no events recorded where the mass is larger 
than roughly ten times its average. We however present now an argument suggesting 
that the tail is faster than exponential, and more particularly that logpi(m) oc — m logm 
when m 3> 1. We first observe that in order to have a large mass in a given cell, 
one needs to transfer to it the mass coming form a large number M of neighboring 
cells. Estimating the probability of having a large mass is equivalent to understand the 
probability of such a transfer. For moving mass from the j'-th cell to the (j — l)-th cell, 
the best configuration is clearly %+i) = (0,1,1). After N time steps with 

this configuration, the fraction of mass transfered is 1 — (1 — r y/2) N . This process is then 
repeated for moving mass to the second neighbor, and so on. After order M iterations, 
the mass in the M-th neighbor is 

1 M 



M=M( m ,N) = 'z i :„ ;: ..rv - < 23 > 



[ 1- (1-7/ 2)"]- 
(1-7/2)' 

This means that 

log [l - m(l - 7/2)*] 
log[l-(l- 7 /2) J 

with the condition that N > — (logm)/[log(l — 7/2)]. The probability of this whole 
process of mass transfer is 

V = [p 2 (l -p)] NM = exp [log(p 2 (l -p))NM(m,N)] . (24) 

All the processes of this type will contribute terms in the right tail of the mass PDF. The 
dominant behavior is given by choosing N = N* such that N* M(m, N*) is minimal. 
Such a minimum cannot be written explicitly. One however notices that, on the one 
hand, if N is much larger than its lower bound (i.e. N ^> — (logm)/[log(l — 7/2)]), 
then N M(m, N) ^> — m (logm)/ [log (1 — 7/2)]. On the other hand when iV is chosen 
of the order of logm, then iV M(m, N) oc m logm. This suggests that the minimum is 
attained for iV* oc logm. Finally, such estimates lead to predict that the right tail of 
the mass probability density function behaves as 

Pi(m) oc exp [— C m logm] , (25) 

where C is a positive constant that depends upon the parameters p and 7. As seen in 
figure El such a behavior is confirmed by numerical experiments. 

The estimations of the left and right tails of the distribution of mass in a given cell 
can be extended to the two-dimensional case. The results do not qualitatively change. 
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Figure 6. Lin-log plot of the one-point PDF of the mass m in one dimension 
represented as a function of m log m for p — 1/2 and various values of the parameter 
7 as labeled; the different colors and symbols are the same as those used in figure |U 
Inset: behavior of the constant C appearing in (|25|) as a function of the ejection rate 
7 for three different values of the fraction of space p occupied by eddies (blue crosses: 
p = 0.1, black times: p = 0.5, red circles: p = 0.9). 



The exponent a of the algebraic behavior at small masses is given as a solution of 
4p3 6p 2 (l — p) 4p(l — p) 2 

(1 _ 7 /4)«+i + (1 - 7 /2) a +! + (1 -3 7 /4) Q+1 + 

+ f^L = ^'P + P 2 )- (26) 

By arguments which are similar to the one- dimensional case and which are not detailed 
here, one obtains also that logpi(m) oc —m logm. Numerical experiments in two 
dimensions confirm these behaviors of the mass probability distribution. As seen from 
figure [7] an algebraic behavior of the left tail of the PDF of m is observed and the value 
of the exponent is in good agreement with f[2"6T) . 



5. Coarse-grained mass distribution 

We investigate in this section the probability distribution of the mass coarse-grained on 
a scale L much larger than the box size £, which is defined as 

i A 

rh L = - 22 m i where K = L / 2t ( 27 ) 

j=-K 

As seen from the numerical results presented on figure El the functional form of the 
PDF pi{fn) is qualitatively similar to that of the mass in a single cell. In particular for 
various values of L it also displays an algebraic tail at small arguments with an exponent 
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Figure 7. Log- log plot of the one-point PDF of the mass in two dimensions for p = 1/2 
and various values of the parameter 7; the symbols and color refer to the same values 
of 7 as in figure HI The integration was done on a domain of 1024 2 cells and time 
average were performed during 3 x 10 4 time steps after a statistical steady state is 
reached. Inset: exponent a of the algebraic left tail as a function of the ejection rate 7 
for three different values of p (blue crosses: p = 0.1, black times: p = 0.5, red circles: 
p = 0.9). The solid lines shows the analytic values obtained from ([26]) . 

which depends both on L and on the parameters of the model. We here present some 
heuristic arguments for the behavior of the exponent. 

For this, we consider the cumulative probability P^(rh) to have mi smaller than 
fh. We first observe that in order to have fh^ small, the mass has to be transfered 
from the bulk of the coarse-grained cell to its boundaries. Assume we start with a mass 
order unity in each of the IK + 1 sub-cells. The best realization to transfer mass is to 
start with ejecting an order-unity fraction of the mass contained in the central cell with 
index j — to its two neighbors. For this the three central cells must have vorticities 
(fi_i, f^o, ^1) = (0,1,0), respectively, during N time steps. After that the second step 
consists in transferring the mass toward the next neighbors; the best realization is then 
to have during N time steps (0_ 2 , fi , tti, f2 2 ) = (0, 1, 1, 1, 0). The transfer toward 
neighbors is repeated recursively. At the j-th step, the best transfer is given by choosing 
Q-j, = (0, 1, 1) and (fij-i, %+i) = (1, 1, 0) during a time N. One can 

easily check that for large N and after repeating K times this procedure, the mass which 
remains in the 2K + 1 cells forming the coarse-grained cell is mi ~ (1 — r y/2) N . The 
total probability of this process is [p 4 (l — p) 2 } ,which leads to estimate the cumulative 
probability of % as 

P<(fh) cx m aL , with a L « \ L log [p\l - pf] . (28) 

2 log(l - 7/2) 

This approach guarantees that the probability density function pi,{fn) = dP^/dfh of 
the coarse-grained mass fh behaves as a power law at small values. Note that only 
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the contribution from realizations with an optimal mass transfer is here evaluated and 
the actual value of the exponent should take into account realizations of the vorticity 
which may lead to a lesser mass transfer. However we expect the estimation given by 
( 128]) to hold for L sufficiently large, because the contribution from realizations with a 
sub-dominant mass transfer become negligible in this limit. 

As to the right tail of pi,(m), one expects a behavior similar to that obtained in 
the case of the one-cell mass distribution, namely logpi(m) oc — rh logm for fh ^> 1. 
Indeed, the probability of having a large mass in a coarse-grained cell should clearly be 
of the same order as the probability of having a large mass in a single cell. This, together 
with the estimates (1281) for the exponent of the left tail, gives a motivation for looking, 
at least in some asymptotics, for possible rescaling behaviors of Phiyn) as a function of 
L and of the ejection rate 7. For instance one can argue whether the limits L — > 00 
and 7 — > are equivalent. The estimation ( 1281) suggests that the exponent c*l depends 
only on the ratio k = L/log(l — 7/2). Note that the limit of small 7 should mimic 
that of small response time of the heavy particles. Rescaling of the distribution of the 
coarse-grained mass was observed in direct numerical simulations of heavy particles in 
turbulent homogeneous isotropic flows [15] . 




Figure 8. Log- log plot of the PDF of the coarse-grained mass mi, in one dimension 
for p = 2/3 and various values of 7 and L associated to three different ratios 
k = Lj log(l — 7/2) as labeled. 

Such a rescaling is confirmed by numerical simulations. Figure M represents the 
PDF of the coarse-grained mass rhi for various values of L and 7 chosen such that the 
ratio k is 8, 10 or 14. While the left-tail and the core of the distribution are clearly 
collapsing, the rescaling is much less evident for the right-hand tail. Obtaining better 
evidence would require much longer statistics in order to resolve the distribution at large 
masses. 
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6. Conclusions 

We introduced here a simple model for the dynamics of heavy inertial particles in 
turbulent flow which solely accounts for their ejection from rotating structures of the 
fluid velocity field. We have shown that this model is able to reproduce qualitatively 
most features of the particle mass distribution which are observed in real turbulent flows. 
Namely the probability density of the mass in cells is shown to behave as a power-law 
at small arguments and to decrease faster than exponentially at large values. Moreover, 
we studied how this distribution depends on the parameters of the model, namely the 
ejection rate of particles from eddies and the fraction of space occupied by them. Such 
dependence reproduce again qualitatively observations from numerical simulations in 
homogeneous turbulent flows. Finally, we have seen that coarse-graining masses on 
scales larger than the cell size is asymptotically equivalent to decrease the ejection rate 
related to particle inertia. This gives evidence that there exists a scaling in the limit of 
large observation scale and small response time of the particles, even if the flow has no 
scale invariance. 

There are several extensions that need to be investigated in order to gain from the 
study of such models a more quantitative information on the distribution of particles 
in real flows. The most significant improvement is to give a spatial structure to the 
fluid velocity. This can be done by introducing a spatial correlation between the 
vorticities of cells. Preliminary investigations suggest that such a modified model could 
be approached by taking its continuum limit. Another effect that may be worth taking 
into consideration is random sweeping of structures by the fluid flow. We assumed that 
the eddies are frozzen (and occupy the same cell) during their whole life time. The 
model could be extended by adding to the dynamics random hops between cells of these 
structures. Another extension could consist in investigating in a systematic manner 
three-dimensional versions of the model. As stated above, many statistical quantities 
may depend on how ejection from rotating regions is implemented in three dimensions. 
Finally, it is worth mentioning again that the main advantage of such models is to give a 
heuristic understanding of the relations between the properties of the fluid velocity field 
and the mass distribution of particles. This first step is necessary in the development of 
a phenomenological framework for describing the spatial distribution of heavy particles 
in turbulent flows. This would allow for using Kolmogorov 1941 dimensional arguments 
to understand how the particle dynamical properties depend on scale. Moreover, such 
a framework could be used to obtain refined predictions accounting for the effect of the 
fluid flow intermittency and to describe the dependence upon the Reynolds number of 
the spatial distribution of particles. 
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